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ABSTRACT 

We present values for angular diameter, flux and extinction for 70 Galactic planetary nebulae 
observed using narrow band filters. Angular diameters are derived using constant emissivity 
shell and photoionization line emission models. The mean of the results from these two mod- 
els are presented as our best estimate. Contour plots of 36 fully resolved objects are included 
and the low intensity contours often reveal an elliptical structure that is not always apparent 
from FWHM measurements. Flux densities are determined, and for both Ha and [O III] there 
is little evidence of any systematic differences between observed and catalogued values. Ob- 
served Ha extinction values are determined using observed Ha and catalogued radio fluxes. 
Ha extinction values are also derived from catalogued Ha and H/3 flux values by means of 
an Ry dependent extinction law. Ry is then calculated in terms of observed extinction values 
and catalogued Ha and H/3 flux values. Comparing observed and catalogue extinction values 
for a subset of Bulge objects, observed values tend to be lower than catalogue values calcu- 
lated with Ry = 3.1. For the same subset we calculate (i?v) = 2.0, confirming that toward 
the Bulge interstellar extinction is steeper than Ry = 3.1. For the inner Galaxy a relation with 
the higher supernova rate is suggested, and that the low-density warm ionized medium is the 
site of the anomalous extinction. Low values of extinction are also derived using dust models 
with a turnover radius of 0.08 microns. 
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1 INTRODUCTION 

The angular diameter and flux of a planetary nebula (PNe) are much 
needed parameters when modelling the optical spectrum and evo- 
lution of the nebula and its central star. From these parameters the 
density and mass of the PNe can be calculated, as well as the ex- 
tinction along the line of sight. PNe also provide an independent 
tracer of stellar populations and extinction in regions where visual 

tracers, such as hot luminous stars, are lacking. 

At the distance of the Bulge of R = 8 kpc iReidl 1993b many 
PNe are only a few arc seconds across. Therefore determination of 
angular diameters, by deconvolving FWHM measurements, need 
to use methods that take into acc ount ionisation structure a nd the 
particular emission line observed. Bedding & Ziilstra 1 1994, here- 
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after BZ) and van Hood (2000, hereafter VH) provide two suitable 
methods of deconvolution. 

Spectroscopic flux measurements can lead to a possible un- 
derestimation if the slit is not properly aligned with the source. 
Therefore accurate determination of fluxes can be achieved using 
calibrated narrow-band filters, provided that correction for the un- 
derlying continuum and any unwanted lines is applied. Extinction 
can then be calculated from observed Ha and catalogue radio flux 
ratios, and compared with values calculated from catalogue Ha/H/3 
flux ratios. A method for determining Ry, the ratio of total to se- 
lective extinction, can also be derived. 

Correction for interstellar extinction is necessary for determin- 
ing emission line ratios, absolute fluxes, and some distance calcu- 
lations. Usually, a single Galactic extinction law is used for dered- 
dening. However, the emission mechanism of ionized media is well 
understood, and PNe are intrinsically very bright. This allows one 
to use PNe to determine extinction constants for Galactic dust dis- 
tribution studies, for line of sights where few other good tracers are 
available. 
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Table 1. Filter characteristics (A). 



Filter 


ESO# 


Peak 


Centre 


FWHM 


jTdX 


Ha on-band* 


654 


6546 


6549 


32 


11 


Ha off-band 


598 


6675 


6673 


65 


37 


[0 III] on-band* 


589 


5016 


5014 


55 


36 


[Om] off-band 


591 


5100 


5110 


60 


43 



"Calculated from observed spectra. 

In order to improve the accuracy of available PNe data, we 
have obtained narrow-band Ha and [O III] CCD images of a sample 
of Galactic PNe and calculated values for angular diameter, flux 
and extinction. 
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2 OBSERVATIONS AND DATA PROCESSING 

We observed 70 Galactic PNe during the two nights of 2002 June 
24 & 25, using the ESO 3.5-m New Technology Telescope (NTT) 
in Chile. We used the EMMI camera (ESO Multi-Mode Instru- 
ment), which at the time had a 2046 x 2046 pixel Tektronix CCD 
with an image scale of / .'27/pixel. For each object a pair of expo- 
sures were taken using on- and off-band Ha filters. 16 of the ob- 
jects were also observed with on- and off-band [O III] filters. The 
FWHM of the on-band filters are adapted to the high velocity range 
of the observed PNe (200 to 300 km s" 1 in the Galactic bulge, cor- 
responding to ~ 5A). The characteristics of the filters, in A, are 
listed in Table Q and plotted in Fig. Q It should be noted that, to 
avoid contamination by the [N II] 6584A line, the Ha on-band fil- 
ter is centred on 6549A and not 6563A. 

In Fig.Qthe dashed curves are the filter response values pro- 
vided by ESO, and the solid curves are response values calculated 
from calibration spectra taken at the NTT (Ha in 2004 February 
using Grism #3 at a resolution of 1 .5A, interpolated to 0. 1 A; [O III] 
in 2003 June using Grism #3 at a resolution of 3A, interpolated to 
lA). The solid vertical lines denote La Silla-air wavelengths and 
the dashed vertical lines the spread of observed wavelengths due to 
Doppler shifting. Assuming the same equivalent widths, the Ha on- 
band filter #654 has blue-shifted by ~ 5A and the [O III] on-band 
filter #589 has red-shifted by ~ 7 A. It can be seen that for [O III] 
this shift makes little difference to the calculation of actual flux 
values. However, for Ha the shift is significant, as actual observed 
wavelengths are in the steep redward side of the curve (especially 
for positive radial velocities relative to the telescope). 

Because the filters in the red arm of EMMI are used in the par- 
allel beam, there is a large change in central wavelength with posi- 
tion on the detector. The solid curve is given for the optical centre of 
the field. Because objects were not always positioned at the optical 
centre of the CCD array (up to ±20 pixels on the first night and up 
to ±50 pixels on the second), wavelength corrections were calcu- 
lated to take into account the filter angle of 7.5°relative to the CCD 
array. This only affected offsets in the j/-axis, with a 20 pixel offset 
generating a wavelength shift of 0.6A. For these reasons, in partic- 
ular, we chose to calculate precise Doppler wavelength shifts when 
determining actual observed flux values. We have not been able to 
determine the reason for the shift in filter transmission curves, so 
would therefore recommend always taking calibration spectra from 
the optical centre of the field in order to check filter response data. 

To reduce the CCD readout time and hence improve 
the observing duty cycle, only the central 400 x 400 pixels 
(108" x 108") of the detector were read out. Calibration of 
the CCD images followed standard procedures: subtraction of a 
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Figure 1. Transmission curves for on-band Ha and [O III] filters. Dashed 
curves are filter response values provided by ESO and solid curves are those 
based on observed calibration spectra. Solid vertical lines indicate La Silla- 
air wavelengths, with dashed vertical lines indicating the spread of observed 
wavelengths due to Doppler shifting. 



constant value to represent the bias, and removal of the back- 
ground continuum, achieved by subtracting the off-band image, 
scaled by the integrated transmission ratio of the two filters 
( J T on dX I J T ff dA), from the on-band image. The ratios were 
0.31 for the Ha filter and 0.84 for the [O III] filter. Where neces- 
sary and prior to this subtraction, pairs of images were brought into 
alignment by integer pixel shifts of the off-band image in the x 
and/or y direction. For the first night's observations this involved 
shifts of only one or two pixels. However, due to taking spectra 
between on- and off-band exposures, the second night's observa- 
tions required shifts of between 20 and 30 pixels. One off-band 
image (356. 1±02.7 Th3-13) had to be rotated as well, to bring it 
into alignment. This procedure removed the background field stars, 
except for objects with spectral lines within the filter band. Cor- 
rection with dome flats was not considered necessary due to the 
sources being very small and pixel to pixel variations in the centre 
of the CCD contributing a possible error in flux values of only 1 
per cent. 

All exposures were for 30 seconds with the autoguider deac- 
tivated, as the NTT can track without guiding for several minutes. 
The resulting images show no sign of elongation of the field stars. 
Ten second exposures of two tertiary standard stars were also taken 
for flux calibration. Seeing conditions varied between one and two 
arc seconds. Subsequent image analysis was done using the ESO 
Munich Image Data Analysis System (MIDAS). 

Bulge PNe were concentrated on, but non-Bulge objects were 
also observed to make best use of telescope time. Using a Galactic 
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Table 2. The parameters Ivan Hool2Q00l tables 1 & 2) for calculating con- 
version factors using equationl2l 



Model 


ai 


a2 


33 


Constant emissivity shell 0.8 


0.3893 


0.7940 


1.2281 


Hcv emission line shell 


0.3605 


0.7822 


1.4937 


[O III] emission line shell 


0.4083 


0.7892 


1.3176 



Bulge criteria for the galactic coordinates (PN G) of | Z < 15° and 
\b\ < 15° gives 56 Bulge PNe and 14 non-Bulge PNe. However, 
adopting the criteria o f lBensbv & Lundstrod J200ll) : |Z| < 10° and 
|6] < 10°; Of! 9 < 9 < 20"; 4.2 mjy < 5 5 gh z < 45 mjy, gives 
39 Bulge PNe and 3 1 non-Bulge PNe. The non-Bulge PNe include 
two Sagittarius dwarf galaxy objects (Dudziak et al. 2000). 



3 DIAMETER DETERMINATIONS 
3.1 Method 

The images presented here are broadened by the seeing. This must 
be taken into account when deriving angular diameters and we now 
describe the methods used to deconvolve the point spread function 
(beam size). 

Using the method of BZ, the point spread function (PSF) was 
approximated with a Gaussian FWHM of 8psf, by fitting Gaussian 
profiles to four distinct field stars in each of the on-band images, 
taking the mean value and expressing the error as the standard de- 
viation of the four readings. The results are listed in column four of 
Table|4| with an overall mean value for PS f of l'.'55. The FWHM 
of the nebulae #fwhm, was measured by fitting two dimensional 
Gaussian profiles to the corrected on-band images. As BZ point 
out, if the objects were Gaussian, deconvolution would be trivial, 
and the true FWHM of the nebulae would simply be 

6» gauss = \A0fwhm) 2 - (0 PS f) 2 . (1) 

Although 6lg auss is not a realistic model for PNe, a correc- 
tion factor 7, can be applied in order to determine the true diam- 
eter, (?true — 7#gauss- Various methods of calculating this fac- 
tor have been investigated by VH and the conversion factor 7 is 
found to be a function of j3, the resolution of the observation, i.e. 
the ratio of the Gaussian diameter and the point spread function, 
ft = # g auss/#PSF- According to VH, 7(/3) can be approximated 
extremely well by the following analytic function (VH method A): 

^ = T4^ +0 " (2) 
which will be used throughout this paper (appropriate parameters 
for equation|2|are given in Table|2j. 

VH finds that the method used by BZ for the calculation of 7 
for disk, sphere and shell geometries is essentially identical to his 
own method B, with results that are also in good agreement with 
his method A (see VH fig. 3), apart from a slight overestimation of 
7 for larger values of j3. In the case of a constant emissivity shell, 
VH fig. 2 shows that the shape of the curve for the conversion fac- 
tor 7 does not change very much as a function of ri/r s , but that 
the height of the curve does. Therefore, the correct value for 7 de- 
pends critically on the assumed value for the inner radius, n, with 
differences of up to 40 per cent, depending on the choice r;. VH 
also examines what he considers to be more realistic geometries 
based on a photoionization model of PNe. This model takes into 
account emissivity and absorption of specific lines as a function of 
distance to the star. As VH point out, the surface brightness profile 



is very sensitive to optical depth effects and thus different emission 
lines will yield different diameters. Measured diameters will also 
be different, depending on which line is used, because of ionisa- 
tion stratification effects. VH also present an alternative algorithm 
for determination of the intrinsic FWHM of a source, using the ob- 
served surface brightness distribution and the FWHM of the beam 
(PSF). This was not investigated as knowledge of the intrinsic sur- 
face brightness distribution is still needed to convert the intrinsic 
FWHM of the source to the true diameter. The reader is referred to 
VH for a full discussion of the various methods that can be used to 
measure angular diameters of extended sources. 

Based on the above we chose to calculate diameters using two 
models: a constant emissivity shell, 0„heii (where an inner diam- 
eter of 0.8# s hcii was chosen because it represents a filling factor 
of 0.5); and photoionization line emission, 6n ue . In Table|4]we list 
deconvolved Gaussian FWHM values, 9 sauss and PSF (beam size) 
values, #psf, so that calculations based on alternative intrinsic sur- 
face brightness distributions can be made. 

3.2 Image Restoration 

Since the seeing was not exceptional during the observations it was 
considered desirable to improve the morphological information of 
the images by restoring them. The critical component of restora- 
tion is fidelity of the PSF, and for that reason only those images 
were selected which had a bright well separated star image on the 
frames. Therefore images with many stars, or with no bright stars, 
were not deemed suitable for restoration. 18 pairs of Ha on- and 
off-band image s were restored u sing the Richardson-Lucy algo- 
rithm I Richardson 1972t lLucvll974l) as implemented in the IRAF 1 
stsdas. analysis. restore package. The images were background sub- 
tracted by a constant and the PSF was cut from the image. Each 
image was restored either for 50 iterations or until convergence 
(X 2 ~ !)■ Using these on-band restored images, diameters for 
#rcstPSF, ^rost gauss, #r CS t shell and 6 rest line were then calculated 
using the same methods as above. The results are listed in Table|3| 
with # rcs t the mean of lest shell and 9 Icst line- Diameters derived di- 
rectly from the on-band-minus-off-band images (Table|4j are also 
included for ease of comparison. An attempt to generate restored 
continuum-free images was not made by subtracting the restored 
off-band image from the restored on-band image. Such difference 
images would likely suffer from high frequency noise, since the re- 
stored PSFs would not be identical in both images and therefore 
could not be used for flux determination. 

3.3 Results 

From the results listed in Table|4|it can be seen that the ionization 
stratification model, 9u nc , consistently yields larger diameter val- 
ues than the constant emissivity shell 0.8 model, (? s helb Calculating 
#sheiio.5 (n = 0.5, filling factor 0.9) yields values within a few per 
cent of 9u nc . However, we believe O s hel\ with a filling factor of 0.5 
(n = 0.8) to be more reasonable. Therefore we present the mean 
value, 9 mcan , of the two models as our best estimate, with #ii ne and 
#shcii the upper and lower limits. For S me an, both axis diameters 
are shown where the axis ratio > 1.1 and the axis difference > 0.3 
arc seconds. Fig.|2|shows that calculated [O III] 9y lnc diameters are 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, Inc., 
under contract to the National Science Foundation of the United States. 
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Table 3. Ha angular diameters derived from 18 restored on-band images, compared with diameters derived directly from on-band-minus-off-band images 
(Tablel4*l. Images restored with the Richardson-Lucy algorithm (Richardson 1972; Lucy 1974). Diameters are in arc seconds and are defined in sectionlSI 



PNG 
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Pvnst pop 

i Lh l r or 


$rcst g 




6 




rest _lu /o 


f?PSF 


e 

gauss 





^10% 


000.3-04.6 


M2-28 


1.04 


4.27 


X 


3.60 


5.9 


X 5.0 


7.5 X 7.0 


1.57 


4.10 


5.8 


9.0 x 8.0 


000.3-02.8 


M3-47 


0.82 


5.63 


X 


3.60 


7.7 


X 5.0 


6.0 x 2.5 


1.10 


5.39 X 3.93 


7.4 X 5.5 


9.0 x 8.0 


000.6-01.3 


B13-15 


0.79 


1.87 


X 


1.55 


2.7 


X 2.3 


5.0 x 4.0 


1.50 


2.27 X 2.03 


3.4 x 3.1 


6.0 x 4.5 


002.3-03.4 


H2-37 


0.72 


3.88 


X 


1.50 


5.3 


X 2.2 


5.5 x 3.0 


1.88 


3.83 X 2.04 


5.5 X 3.2 


7.0 x 5.5 


002.8-02.2 


Pe2-12 


0.85 


6.73 


X 


1.89 


9.2 


X 2.7 


7.0 X 3.0 


2.00 


7.01 X 2.14 


9.8 x 3.3 


10.5 x 5.0 


003.0-02.6 


KFL4 


0.63 


1.19 


X 


1.04 


1.7 


x 1.5 


2.5 


1.81 


1.46 


2.2 




003.7-04.6 


M2-30 


0.68 


2.81 


X 


2.81 


3.9 


x 3.9 


4.0 


1.57 


3.00 X 2.58 


4.4 x 3.8 




004.8-05.0 


M3-26 


0.70 


6.11 


X 


5.84 


8.4 


x 8.0 


10.0 X 8.5 


1.54 


5.93 


7.7 


11.0 x 9.5 


005.2-18.6 


StWr2-21 


0.85 


0.74 


X 


0.76 


1.2 


x 1.2 


2.5 X 2.0 


1.80 


1.58 x 1.21 


2.5 X 2.0 




008.6-07.0 


He2-406 


1.16 


3.73 


X 


2.57 


5.2 


x 3.7 


6.5 X 5.5 


1.56 


4.06 X 3.37 


5.8 x 4.9 


8.0 x 7.5 


283.8+02.2 


My60 


0.80 


6.04 


X 


5.60 


8.3 


X 7.7 


10.0 


1.65 


5.65 


7.9 


11.0 X 10.5 


346.3-06.8 


Fg2 


0.46 


3.39 


X 


2.93 


4.7 


x 4.0 


5.5 x 4.5 


1.17 


3.44 x 2.96 


4.9 x 4.2 


6.5 x 5.5 


352.6+00.1 


Hl-12 


1.18 


6.29 


X 


5.19 


8.7 


X 7.2 


10.0 x 8.5 


1.13 


5.76 


7.9 


12.0 x 10.0 


352.6+03.0 


Hl-8 


0.51 


1.80 


X 


1.50 


2.5 


x 2.1 


4.0 x 3.0 


1.36 


2.10 x 1.73 


3.1 x 2.6 




352.8-00.2 


Hl-13 


0.57 


9.88 


X 


9.03 


13.5 


x 12.3 


13.5 x 12.0 


1.57 


9.27 


12.7 


13.5 X 12.0 


354.9+03.5 


Th3-6 


1.28 


1.74 


X 


1.78 


2.6 


x 2.7 


3.0 


2.38 


2.64 X 2.34 


4.1 x 3.7 




357.1-06.1 


M3-50 


0.49 


4.81 


X 


0.62 


6.6 


x 0.9 


7.0 X 2.0 


1.25 


4.82 X 1.07 


6.7 x 1.7 


8.5 X 3.5 


357.1-04.7 


Hl-43 


0.72 


0.54 


X 


0.34 


0.9 


x 0.6 


2.0 


1.71 


1.05 X 0.75 


1.7 x 1.3 





in most cases smaller than those for Ha, as would be expected. As 
Table|4]shows, 36 of the objects are fully resolved (#ii nc > 2#psf), 
and Fig.[3]plots 26" x 26" images of these objects, with contours 
expressed as percentages of the peak nebula flux. The morphology 
of the nebulae are apparent and the structures seen in the images 
(elliptical, disc, bipolar, etc.) c an be compared with the evolution- 
ary sequences of lBalicldil98% . 

25 objects were sufficiently resolved so that a direct diameter 
measurement could be made. For PNe with a well defined outer 
radius (i.e., with a steep drop-off to zero surface brightness at a cer- 
tain radius) it is easy to measure directly the radius for a specific 
surface brightness. Published data and HST images show that flux 
usually dro ps very fast between c ontour levels in the interval 15 to 
5 per cent (Bensbv & Lundstrom 2001). Therefore 10 per cent of 
the peak surface brightness contour diameter measurements were 
made with an estimated error of ±0'.'5, and these are listed in col- 
umn nine of Table |4] It is worth noting that in these cases the low 
intensity contours seen in Fig.|5|often reveal a larger elliptical struc- 
ture than that apparent from the 2D FWHM axis ratio measurement 
of the bright nebula core (column twelve of Table|4j- For compari- 
son columns ten and eleven of Table|4]list the catalogued radio and 
optical diameters from the compilat i on avai lable in the Strasbourg 
- ESO Catalogue lAcker et alll992t) . lTvlenda et alJiiooH) have re- 
cently redetermined a large set of PNe diameters: we find good 
agreement with their results for all but the most compact (< 2") 
objects. 

Table [5] lists Ha angular diameters derived from the 18 re- 
stored on-band images, together with diameters derived directly 
from the on-band-minus-off-band images (Table |4j- For 12 ob- 
jects column eight ((9 gauss ) also shows both axis values, which are 
only averaged in Table [4] f3 rC st is presented as the mean value of 
#rcst shdi and 6» rcst Unc . Fig. compares angular diameters # mcan 
and # r cst, and it can be seen that, regardless of PNe axis, re- 
stored diameters tend to be around one arc second smaller below 
4 arc seconds. The data suggests the simple linear relationship 

"rest — 

1.10 

mean 

1. Direct diameter measurements, # re stio%> 
were made at 10 per cent of the peak surface brightness contour for 
13 of the restored on-band images, and as Table|5|shows, these are 
consistently lower by 1-2 arc seconds than the non-restored 10 per 
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[Offl] Angular Diameter 6| ire (arc sees) 
Figure 2. Comparison of Ha and [O III] angular diameters £>i; nc . 

cent measurements, f5 10 % (except for 352.8—00.2, Hl-13). It could 
be that the restoration algorithm removes the detail that can be seen 
in the very low intensity contours of some PNe in Fig. [5] 

3.4 Discussion 

By their nature PNe at a given distance are likely to exhibit a range 
of angular diameters depending on where they are in their evolu- 
tionary sequence of returning shed material from the proginator 
star to the ISM. This is in addition to observed angular size dif- 
ferences for the same source, e.g. optical versus radio or specific 
line emission. 

It is apparent from a visual inspection of our images that not 
all PNe have a well-defined outer radius. For such soft boundary 
PNe a 10 per cent radius does not have a direct physical mean- 
ing and usually does not correlate well with converted deconvolved 
FWHM diameters. This is not surprising as 10 per cent (or lower 
if signal-to-noise will allow) radius measurements 'see' into the 
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Figure 3. Ha images with contour levels at 80, 65, 50, 35, 20, 10, (5, 2, 1, 0.5, 0.2, 0.1) per cent of the nebula peak (lowest indicated after PNe name). All 
images are 26" X 26" with north up and east to the left. To filter the noise, the images have been smoothed with a 3 pixel (0^81) box. 




Figure 3. (continued) Ha images with contour levels at 80, 65, 50, 35, 20, 10, (5, 2, 1, 0.5, 0.2, 0.1) per cent of the nebula peak. Actual lowest contour level 
plotted dependent on signal-to-noise ratio, and indicated after PNe name. The images have been smoothed with a 3 pixel (0'.'81) box. 
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Table 4. Angular diameters of 69 planetary nebulae. Diameters are in arc seconds and are defined in sectionl3l 
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Ha 


1.10 + 0.09 


4.67 


5.8 


7.1 


7.4 X 5.5 


±0.6 


9.0 x 8.0 




St. 


1.35 


fr, irreg 


000.6-01.3 


B13-15 


Ha 


1.50 + 0.03 


2.15 


3.0 


3.5 


3.4 X 3.1 


±0.3 


6.0 x 4.5 




3.0 


1.08 


fr, ellip 


000.8-01.5 


BIO 


Ha 


1.30 + 0.02 


1.89 


2.6 


3.1 


2.8 


±0.2 






St. 


1.03 


fr 


002.3-03.4 


H2-37 


Ha 


1.88 + 0.02 


2.98 


4.0 


4.8 


5.5 X 3.2 


±0.4 


7.0 x 5.5 




4.2 


1.54 


fr, ellip 


002.8-02.2 


Pe2-12 


Ha 


2.00 + 0.06 


4.70 


6.1 


6.6 


9.8 X 3.3 


±0.2 


10.5 X 5.0 




5.0 


2.49 


fr, disc 


003.0-02.6 


KFL4 


Ha 


1.81 + 0.04 


1.46 


2.2 


2.3 


2.2 


±0.1 






3.0 


1.01 


pr 


003.6+03.1 


M2-14 


Ha 


1.94 + 0.08 


1.68 


2.5 


2.9 


3.1 x 2.2 


±0.2 




2.2 


St. 


1.17 


pr 


003.7-04.6 


M2-30 


Ha 


1.57 + 0.06 


2.79 


3.7 


4.5 


4.4 X 3.8 


±0.4 




3.5 


9.0 


1.12 


fr 


003.8-17.1 


Hb8 


Ha 


1.50 + 0.09 


0.66 


1.0 


1.2 


1.1 


±0.1 




0.8 


5.0 


1.04 


br 


004.0-11.1 


M3-29 


Ha 


1.64 + 0.07 


5.24 


6.7 


8.0 


7.3 


±0.7 


10.5 x 9.5 


13.0 


8.2 


1.06 


fr, ellip 


004.7-11.8 


He2-418 


Ha 


1.18 + 0.05 


7.40 


9.2 


11.1 


12.8 x 7.6 


± 1.0 


14.0 x 8.5 




13.0 


1.68 


fr, irreg 


004.8-22.7 


He2-436 


Ha 


1.38 + 0.14 


0.32 


0.5 


0.6 


0.6 


±0.0 






< 10 


1.05 


br 


004.8-05.0 


M3-26 


Ha 


1.54 + 0.02 


5.93 


7.5 


8.0 


7.7 


±0.3 


11.0 X 9.5 




8.6 


1.03 


fr, disc 


005.2-18.6 


StWr2-21 


Ha 


1.80 + 0.13 


1.40 


2.1 


2.4 


2.5 X 2.0 


±0.2 






< 5 


1.10 


pr 


006.0-41.9 


PRMG1 


Ha 


1.38 a 


2.65 


3.5 


4.2 


4.4 X 3.3 


±0.3 






8.2 


1.29 


fr, irreg 


006.8-19.8 


Wrayl6-423 


Ho 


1.78 + 0.12 


1.08 


1.6 


1.9 


1.9 x 1.6 


±0.1 






St. 


1.05 


P r 


008.6-07.0 


He2-406 


Ha 


1.56 + 0.04 


3.72 


4.8 


5.8 


5.8 x 4.9 


±0.5 


8.0 X 7.5 




3.0 


1.17 


fr, irreg 


008.6-02.6 


MaCl-11 


Ha 


1.82 + 0.04 


1.97 


2.8 


3.3 


3.3 x 2.9 


±0.2 






St. 


1.09 


P r 


009.4-09.8 


M3-32 


Ha 


1.52 + 0.04 


4.32 


5.5 


6.7 


6.8 x 5.4 


±0.6 


8.5 x 7.5 


7.5 


6.0 


1.24 


fr, ring 


037.5-05.1 


A58 


[Oin] 


1.77 + 0.05 


0.70 


1.1 


1.2 


1.4 x 0.9 


±0.0 








1.06 


br, irreg 


283.8+02.2 


My60 


Ha 


1.65 + 0.02 


5.65 


7.2 


8.6 


7.9 


±0.7 


11.0 x 10.5 




7.6 


1.07 


fr 






[Om] 


1.97 + 0.03 


6.43 


8.2 


8.8 


8.5 


±0.3 


12.5 x 12.5 






1.07 


fr 


289.8+07.7 


He2-63 


Ha 


1.47 + 0.06 


1.83 


2.6 


3.0 


3.0 x 2.6 


±0.2 






3.0 


1.11 


fr 






[Om] 


1.55 + 0.03 


1.87 


2.6 


2.8 


3.1 x 2.6 


±0.1 








1.12 


P r 


292.8+01.1 


He2-67 


Ha 


1.55 + 0.01 


2.18 


3.0 


3.6 


3.6 x 3.0 


±0.3 


6.0 x 4.5 




< 5 


1.16 


fr, ellip 


296.4-06.9 


He2-71 


[Om] 


2.04 + 0.20 


0.36 


0.6 


0.6 


0.6 


±0.0 








1.08 


br 


305.1+01.4 


He2-90 


Ha 


1.39 + 0.04 


0.48 


0.8 


0.9 


0.8 


±0.1 






< 10 


1.01 


br 






[Om] 


1.71 + 0.01 


0.63 


1.0 


1.1 


1.3 x 0.8 


±0.0 








1.06 


br 


315.4+09.4 


He2-104 


Ha 


1.73 + 0.05 


0.40 


0.6 


0.7 


0.7 


±0.0 






< 5 


1.02 


br, outflow 


321.3-16.7 


He2-185 


Ha 


1.73 + 0.06 


1.53 


2.2 


2.6 


2.8 X 2.0 


±0.2 






< 5 


1.17 


pr 


322.4-00.1 


Pe2-8 


Ha 


1.68 + 0.06 


0.97 


1.5 


1.7 


1.9 x 1.3 


±0.1 






1.6 


1.09 


pr 


345.0-04.9 


Cnl-3 


Ha 


1.48 + 0.02 


0.74 


1.2 


1.3 


1.2 


±0.1 






< 5 


1.02 


br 






[Om] 


1.75 + 0.01 


0.71 


1.1 


1.2 


1.2 


±0.0 








1.02 


br 


345.0+03.4 


Vdl-4 


Ha 


1.79 + 0.03 


0.66 


1.0 


1.2 


1.1 


±0.1 






5.0 


1.03 


br 






[Om] 


1.62 + 0.03 


0.61 


1.0 


1.0 


1.0 


±0.0 








1.03 


br 


345.0+04.3 


Vdl-2 


Ha 


1.86 + 0.01 


0.76 


1.2 


1.4 


1.3 


±0.1 






St. 


1.00 


br, diffuse 






[Om] 


1.70 + 0.03 


0.79 


1.2 


1.3 


1.3 


±0.0 








1.70 


br 


346.0+08.5 


He2-171 


Ha 


1.01 + 0.03 


0.08 


0.1 


0.2 


0.1 


±0.0 






< 10 


1.04 


br 






[Om] 


1.10 + 0.07 


0.28 


0.5 


0.5 


0.5 


±0.0 








1.00 


br 


346.3-06.8 


Fg2 


Ha 


1.17 + 0.04 


3.20 


4.1 


5.0 


4.9 X 4.2 


±0.4 


6.5 x 5.5 




< 5 


1.14 


fr 






[Om] 


1.13 + 0.03 


3.17 


4.1 


4.4 


4.8 X 4.2 


±0.1 








1.13 


fr 


347.4+05.8 


Hl-2 


Ha 


1.14 + 0.02 


0.55 


0.9 


1.0 


1.1 X 0.8 


±0.1 




0.8 


St. 


1.07 


br 






[Om] 


1.26 + 0.03 


0.59 


0.9 


1.0 


1.2 X 0.8 


±0.0 








1.07 


br 


348.0+06.3 


MGP1 


Ha 


1.98±0.14 c 


4.92 


6.4 


7.7 


7.6 x 6.4 


±0.6 


12.0 x 9.5 


3.8 




1.17 


fr, ellip 


348.8-09.0 


He2-306 


Ha 


1.40 + 0.02 


2.04 


2.8 


3.3 


3.1 


±0.3 






3.0 


1.07 


fr 






[Om] 


1.57 ± 0.01 


2.01 


2.8 


3.0 


2.9 


± 0.1 








1.06 


pr 


349.8+04.4 


M2-4 


Ha 


1.12 + 0.02 


1.40 


2.0 


2.3 


2.1 


±0.2 




2.0 


< 5 


1.09 


fr, ellip 






[Om] 


1.21 + 0.03 


1.32 


1.9 


2.0 


2.0 


±0.1 








1.07 


Pr 


350.8-02.4 


Hl-22 


Ha 


2.13 + 0.09 


2.17 


3.1 


3.7 


3.4 


±0.3 






St. 


1.01 


pr 


350.9+04.4 


H2-1 


Ha 


1.38 + 0.03 


1.79 


2.5 


2.9 


2.7 


±0.2 




2.2 


5.6 


1.00 


fr 






[Om] 


1.32 + 0.04 


0.78 


1.2 


1.3 


1.2 


±0.0 








1.01 


br 


351.1+04.8 


Ml-19 


Ha 


1.32 + 0.02 


1.99 


2.7 


3.2 


3.2 X 2.7 


±0.3 




2.6 


8.0 


1.13 


fr 


351.2+05.2 


M2-5 


Ha 


1.33 + 0.01 


3.41 


4.4 


5.3 


4.9 


±0.4 


6.5 x 6.5 


5.0 


5.0 


1.03 


fr, outflow 


351.3+07.6 


Hl-4 


Ha 


1.51 + 0.02 


0.64 


1.0 


1.2 


1.1 


±0.1 






St. 


1.02 


br 


351.9-01.9 


Wray 16-286 


Ha 


1.69 + 0.06 


1.16 


1.8 


2.0 


2.2 x 1.6 


±0.1 






St. 


1.12 


pr, ellip 


351.9+09.0 


PC13 


Ha 


2.07 + 0.04 


5.15 


6.7 


8.0 


8.6 x 6.0 


±0.7 


10.0 x 8.5 




7.0 


1.39 


fr, disc 


352.0-04.6 


Hl-30 


Ha 


1.78 + 0.05 


1.63 


2.4 


2.8 


3.0 x 2.1 


±0.2 






5.4 


1.20 


pr, ellip 



1 Error expressed as the standard deviation of four 0psf values except: "one value, 'two values, c three values. 

2 Maximum/minimum axis ratio of planetary nebulae from 2D FWHM measurements. 

3 br, barely resolved, 9\i nc < 9pgp; pr, partly resolved, 9psF < #iinc < 26>pgp; fr, fully resolved, 6>i; nc 26>pgp. 
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Table 4. (continued) Angular diameters of 69 planetary nebulae. Diameters are in arc seconds and are defined in sectionl3l 



PNG 


Object 


Line 




SpSF 1 


Q 

gauss 




^linc 







1U /o 


^rad 


Qopt 


AR 2 


Notes 3 


352.1+05.1 


M2-8 


Ha 


1 


.10 + 0.01 


2.55 


3.3 


4.0 


3.8 X 3.5 


±0.3 


5.0 x 5.0 


3.7 


4.2 


1.10 


fr, outflow 






[OlII] 


1. 


.24 + 0.04 


2.37 


3.2 


3.4 


3.3 


±0.1 








1.05 


fr, outflow 


352.6+00.1 


Hl-12 


Ha 


1. 


.13 + 0.03 


5.76 


7.2 


8.7 


7.9 


±0.8 


12.0 x 10.0 


11.0 


6.8 


1.00 


fr 


352.6+03.0 


Hl-8 


Ho 


1 


.36 + 0.02 


1.92 


2.6 


3.1 


3.1 x 2.6 


±0.2 






3.4 


1.14 


fr 


352.8-00.2 


Hl-13 


Ha 


1. 


.57±0.16 c 


9.27 


11.5 


14.0 


12.7 


± 1.2 


13.5 x 12.0 


14.0 


9.6 


1.08 


fr, ring 


354.2+04.3 


M2-10 


Ha 


1. 


.64 + 0.03 


2.93 


3.9 


4.7 


5.2 X 3.4 


±0.4 


6.5 x 5.5 


4.0 


4.0 


1.42 


fr, irreg 


354.5+03.3 


Th3-4 


Ha 


1 


.84 + 0.11 


0.97 


1.5 


1.7 


2.0 X 1.2 


±0.1 






St. 


1.12 


br 


354.9+03.5 


Th3-6 


Ha 


2 


.38 + 0.08 


2.49 


3.6 


4.2 


4.1 X 3.7 


±0.3 






St. 


1.07 


pr 


355.1-06.9 


M3-21 


Ha 


1. 


.51 + 0.02 


1.00 


1.5 


1.8 


1.8 X 1.5 


±0.1 






< 5 


1.08 


pr 


355.1+02.3 


Th3-ll 


Ha 


1 


.42 ± 0.04 


1.21 


1.8 


2.1 


2.1 X 1.7 


±0.2 






St. 


1.10 


P r 


355.4-02.4 


M3-14 


Ha 


1 


.52 + 0.02 


2.40 


3.3 


3.9 


3.6 


±0.3 


8.0 x 5.0 


2.8 


7.2 


1.01 


fr, ellip 


355.7-03.5 


Hl-35 


Ha 


1 


.40 ± 0.04 


0.90 


1.4 


1.6 


1.5 


±0.1 




1.1 


2.0 


1.07 


pr 


355.9+02.7 


Th3-10 


Ha 


1 


.24 + 0.03 


1.58 


2.2 


2.6 


2.4 


±0.2 




2.0 


St. 


1.07 


fr 


356.1+02.7 


Th3-13 


Ha 


1 


.63 + 0.02 


0.45 


0.7 


0.8 


1.1 X 0.2 


±0.1 




< 2 


St. 


1.07 


br 


356.8+03.3 


Th3-12 


Ha 


1 


.75 + 0.13 


1.22 


1.8 


2.1 


2.3 x 1.6 


±0.1 




1.2 


St. 


1.12 


pr 


357.1-06.1 


M3-50 


Ha 


1 


.25 + 0.04 


3.06 


4.0 


4.8 


6.7 x 1.7 


±0.4 


8.5 x 3.5 




4.0 


3.03 


fr, disc 


357.1-04.7 


HI -43 


Ha 


1 


.71 + 0.03 


0.91 


1.4 


1.6 


1.7 x 1.3 


±0.1 




3.0 


2.0 


1.08 


br 


357.1+01.9 


Th3-24 


Ha 


1 


.24 + 0.04 


4.61 


5.8 


7.0 


6.4 


±0.6 


8.0 x 7.0 




St. 


1.07 


fr, irreg 


357.1+03.6 


M3-7 


Ha 


1 


.39 + 0.02 


3.15 


4.1 


4.9 


4.9 x 4.1 


±0.4 


6.5 x 5.5 


4.7 


5.8 


1.18 


fr 


357.6-03.3 


H2-29 


Ha 


1 


.45 + 0.02 


7.31 


9.1 


11.0 


12.5 X 7.7 


± 1.0 






4.8 


1.61 


fr, irreg 


358.5+03.7 


A12-B 


Ha 


1 


.56 + 0.04 


1.62 


2.3 


2.7 


3.0 x 2.1 


±0.2 






St. 


1.22 


pr, irreg 


358.5-01.7 


JaSt64 


Ha 


1 


.19 + 0.05 


0.98 


1.5 


1.7 


1.6 


±0.1 








1.04 


pr 


358.7+05.2 


M3-40 


Ha 


1 


.49±0.12 b 


1.96 


2.7 


3.2 


3.0 


±0.3 




2.5 


St. 


1.05 


fr 


359.2+04.7 


Th3-14 


Ha 


1 


.48 + 0.05 


1.15 


1.7 


2.0 


1.9 


±0.1 




< 3.3 




1.03 


pr 



1 Error expressed as the standard deviation of four 8p$F values except: "one value, 'two values, c three values. 

2 Maximum/minimum axis ratio of planetary nebulae from 2D FWHM measurements. 

3 br, barely resolved, 0\i nc < $psf; pr, partly resolved, 0psf < #iine < 2£?psf; fr, fully resolved, 6>i; nc 26>pgp. 



Table 5. Flux and extinction values of 70 planetary nebulae. Observed (equation^ and catalogue (equation llOl extinction are defined in sectionl5l 



PNG 


Object 


Line 


Radial Vel. fkms" 1 ) 


Flux (log mWm" 2 ) 


6 cm 


Ha Extinction 


Rv 


ARy 


Notes 5 








Observed 


Catalogue 


Observed 


Catalogue 


(mJy) 4 


Observed 


Catalogue 








000.3-04.6 


M2-28 


Ha 


-43.5 


-29.9 


-11.70 


-11.73 


10 


0.66 


0.98 


1.99 


0.62 


B 


000.3-02.8 


M3-47 


Ha 


-17.9 


-16.0 


-12.47 


-12.88 






1.58 








000.6-01.3 


B13-15 


Ha 


-114.8 




-12.59 


-12.92 






2.89 








000.8-01.5 


BIO 


Ha 


79.9 




-12.18 


-11.82 






2.04 








002.3-03.4 


H2-37 


Ha 


-129.8 


-156.8 


-11.86 


-12.06 






0.69 








002.8-02.2 


Pe2-12 


Ha 


108.8 




-11.86 


-12.46 


2 


0.13 


1.06 


0.57 


0.41 


lrf 


003.0-02.6 


KFL4 


Ha 


-12.9 


17.0 


-12.59 


-12.81 


1 


0.59 


1.02 


1.71 


0.56 


lrf 


003.6+03.1 


M2-14 


Ha 


-25.0 


-48.2 


-11.35 


-12.15 


39 


0.90 


1.41 


1.87 


0.42 


B 


003.7-04.6 


M2-30 


Ha 


130.0 


154.9 


-11.29 


-11.18 


14 


0.40 


0.74 


1.59 


0.75 


B 


003.8-17.1 


Hb8 


Ha 


-170.1 


-180.5 


-11.30 


-11.28 


6 


0.03 


0.20 


0.64 


2.30 


lrf 


004.0-11.1 


M3-29 


Ha 


26.2 


49.3 


-11.20 


-11.10 


12 


0.25 


0.12 


10.81 


31.50 


B 


004.7-11.8 


He2-418 


Ha 


87.7 


109.3 


-12.05 


-11.90 






0.12 








004.8-22.7 


He2-436 


Ha 


128.4 


133.0 


-11.86 


-11.50 


3.9 


0.41 


0.50 


2.45 


1.36 


lrf 


004.8-05.0 


M3-26 


Ha 


-24.7 


-9.8 


-11.41 


-11.26 


8 


0.28 


0.43 


1.90 


1.40 


lrf 


005.2-18.6 


StWr2-21 


Ha 


136.3 


133.0 


-12.26 


-12.46 






0.32 








006.0-41.9 


PRMG1 


Ha 


-8.1 




-12.95 


-12.77 






0.18 








006.8-19.8 


Wray 16-423 


Ha 


139.2 


133.0 


-11.78 


-11.48 


4.7 


0.42 


0.36 


3.78 


2.47 


lrf 


008.6-07.0 


He2-406 


Ha 


-8.5 


28.2 


-12.12 


-12.38 






0.84 








008.6-02.6 


MaCl-11 


Ha 


-30.8 


-89.0 


-12.33 


-12.25 






1.76 








009.4-09.8 


M3-32 


Ha 


40.1 


58.4 


-11.39 


-11.22 


12 


0.43 


0.53 


2.46 


1.29 


B 


009.8-07.5 


GJJC1 


[Om] 




-32.0 


-12.74 
















037.5-05.1 


A58 


Ha 
[Om] 




70.0 


-13.40 
-12.91 
















283.8+02.2 


My60 


Ha 
[Om] 




> 


-11.53 
-10.68 


-11.04 
-10.64 


60 


< 1.27 


0.67 






npc 



4 Bold values are for PNe in the Sagittarius dwarf galaxy iDudziak etall2 000). All other values from the Strasbourg - ESO Catalogue <Ackeretalll992l) . 

5 B, Bulge subset B; npc, non-photometric conditions; hbt, high brightness temperature, Tb > 1000 K; lrf, low radio flux, Si/ < 10 mJy. 



10 P. M. E. Ruffle et al. 



Table 5. (continued) Flux and extinction values of 70 planetary nebulae. Observed (equationQ and catalogue (equation llOl extinction are defined in sectionl5l 



PNG 


Object 


Line 


Radial Vel. (kms^ 1 


) 


Flux (log mWm 2 ) 


6 cm 


Ho Extinction 


By 


AR V 


Notes 5 








Observed 


Catalogue 


Observed 


Catalogue 


(mjy) 4 


Observed 


Catalogue 








289.8+07.7 


He2-63 


Ha 
[Oin] 




123.2 




-11.92 

11 1 A 


-11.91 

1 1 AQ 


12 


0.96 


0.27 






npc 


292.8+01.1 


He2-67 


Ha 




59.5 




-12.47 


-11.17 


41 


< 1.51 


0.83 






npc 


296.4-06.9 


He2-71 


[Om] 








1 1 71 
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5 B, Bulge subset B; npc, non-photometric conditions; hbt, high brightness temperature, T^, > 1000 K; lrf, low radio flux, S u < 10 mJy. 
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Figure 4. Comparison of Ha angular diameters mea n and 9 TCB 
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extended low emission regions surrounding PNe, whereas FWHM 
determinations tend to reveal the compact core. These differing val- 
ues actually give us alternative views of a PNe's morphology and 
tell us more about the object than a single arbitrary radius value. 



4 FLUX DETERMINATIONS 
4.1 Method 

The total flux detected is reduced by the transmission of the Ha 
and [Om] filters. This and other factors must be taken into ac- 
count when determining actual flux values, and we now describe 
the methods used to arrive at the values listed in Table|5| 

The central position of each object was computed by fitting 
a Gaussian to the image distribution and then the magnitude was 
calculated by integrating over successively larger diameters of the 
central area. By plotting these radii against magnitude it was possi- 
ble to determine the radius at which the flux level matched that of 
the background continuum, and therefore the total flux of the ob- 
ject in ADU counts per 30 second exposure. These profiles were 
compared with those of standard stars LTT1020 and LTT624 to de- 
termine the aperture correction. 

To determine detector flux values in mWm~ 2 (erg cm _2 s _1 ), 
the ratio of ADU count to flux conversion factor was calculated 
by measuring the count in ADUs per ten seconds of LTT1020 
and LTT6248, and comparin g it with tabul ated fluxes for these 
spectrophotometric standards l H amuv et alll994h . The calibration 
counts were factored for the shorter exposure time, and the aperture 
correction of —5 per cent was applied. Apart from a brief period of 
non-photometric conditions caused by cloud (indicated as 'npc' in 
Table observations were considered photometric as the com- 
puted ADU count to flux ratios varied by no more than 2.1 per cent 
between the two night's observations. 

In order to determine actual Ha 6565A flux values, any con- 
tribution from [Nil] lines at 6548A and 6584A (vacuum), passing 
through ESO filter #654, was taken into account by assuming the 
following expression for the total detected flux: 



N 



served/Catalogue HaAFlux Distribution 



0.0-0.1 0.1-0.2 0.2-0.3 0.3-0.4 0.4-0.5 0.5-0.6 0.6-0.7 0.7-0.8 0.8-0.9 0.9-1.0 
AFlux (log mW rrr 2 ) 

Figure 5. Distribution of A_F(Hq), the absolute difference between ob- 
served and catalogued Ha flux values (log mWm — 2 ). 



where FX and TX are actual flux values and filter transmission co- 
efficients (corrected for the appropriate Doppler shift) respectively. 
Using catalogue [N II] to Ha flux ratios i?6584 = ^6584/^6565 an d 
J?6548 = ^6548 / Feses (which avoid a dependence on absolute val- 
ues), the actual Ha 6565A flux can be calculated using: 



F(Ha) 



-Pdet/(T6565 + #658426584 + ReS&TeB 



(4) 



2*656526565 + ^6584 26584 + i*6548 76548 



(3) 



[N II] 65 48 flux values were calculated from the [N II] 6584 catalogue 
values I A cker et alJI 199211 using a factor of 2.97. In the case of 
[O III], the flux observed through the filter was factored for the ap- 
propriate transmission coefficient of the [O III] wavelength, so as to 
arrive at actual observed [O III] flux values. 

The Doppler shifted Ha and [O III] wavelengths observed 
at the telescope were computed from catalogue radial velocities 
iDurand et alJll998t lAcker et alJll992l) less a barycentric correc- 
tion appropriate to the date and time at La Silla, or in the absence 
of a catalogue value, derived from spectra taken during our obser- 
vations. Columns four and five of Table [5] compare catalogue ra- 
dial velocities with observed values, which we consider accurate 
to ±25 km s _1 . For three objects it was not possible to determine 
the radial velocity at all. Using the corrected filter response curves 
(measured at La Silla air pressure), the appropriate transmission 
coefficient was then applied to arrive at actual flux values. La Silla- 
air wavelengths were 6563. 3A for Ha and 5007. 3A for [O III], with 
6548.5A and 6583.9A assumed for [N II] using the same AA as that 
of Ha. Although stellar velocities in the Bulge are typically of the 
order of ±100 km s -1 iMinniti 1996), or about ±2A, we have cho- 
sen to calculate precise Doppler wavelength shifts, in order to take 
into account the sensitivity of the Ha filter transmission curve to 
small changes in observed wavelength. 



4.2 Results and error factors 

Table|5|lists the newly determined Ha and [Om] flux densities in 
column six. For comparison columns seven and eight list Ha and 
[O III] flux values (calculated from catalogued H/3 flux and line 
intensities) and t he 6 cm ra dio flux values from the Strasbourg - 
ESO Catalogue jAcker et alj|l992h . Radio values in bold are for 
the tw o PNe in the Sagittarius dwarf galaxy (from lDudziak et ail 
2000). Fig.|5|shows the distribution of AF(Ha), the absolute dif- 
ference between observed and catalogued Ha flux values. 7 out of 
60 objects are shown to have AF(Ha) > 0.4 (log mWm -2 ). A 
comparison of AF(Ha) with r9 linc shows no evidence of a correla- 
tion between these flux differences and angular diameter. 
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Figure 6. Comparison of observed and catalogued Ha and [O III] flux val- 
ues (log mWm - 2 ). Diamonds denote objects observed on the first night, 
triangles the second. 



Fig.[6]compares observed and catalogued Ha and [O III] flux 
values. Four (non-Bulge) objects (indicated in Table |5} are ex- 
cluded from Figs.|5|&|6|because atmospheric conditions for those 
exposures were non-photometric. It can be seen that there is little 
evidence of any systematic differences between observed and cata- 
logued flux values. 

Possible sources of error in our flux determinations are: poor 
photometry; residual stellar background flux after subtraction of 
the continuum; an off-centred peak in the Ha on-band filter re- 
sponse curve; estimation of the flux contribution from [N II] lines; 
the accuracy of transmission coefficients in the corrected filter re- 
sponse tables (particularly for Ha); change in passbands of filters 
due to their age, input beam and ambient temperature; integrating 
the transmission profile of the filter response curves; uncertainties 
in radial velocities; changes in air pressure (airmass); and the aper- 
ture correction factor. Of these, we consider the accuracy of the 
filter response data the most significant and therefore suggest an 
error factor in the calculated flux values (mWm -2 ) of ±5 per cent 
for blue-shifted objects and ±10 per cent for red-shifted objects. 



5 EXTINCTION DETERMINATIONS 
5.1 Method 

Observed Ha flux values are reduced from emitted flux values by 
interstellar extinction. By comparing these observed values with 
catalogued optical and radio fluxes, extinction values can be de- 
termined, and we now describe the methods used to arrive at the 
values list ed in TableEl 

From lPottaschl (l984, hereafter SP) the radio to H/3 flux ratio 



is S„/F(Kp) = 2.51 x 10 7 T e u 



Y (Jy/mWrn -2 ) where 



T c is the electron temperature in K, v is the radio frequency in 
GHz and Y is a factor incorporating the ionized He/H ratio. Setting 
T c = 10 4 K, v = 5 GHz, Y — 1.1 and converting to mJy gives 



S„(6cm)/F(H/3) = 3.10 x 10 12 (mJy/mWm" 2 ). 



(5) 



H/3 extinction (Cp) can be defined as the ratio of expected to ob- 
served H/3 flux: 



Cp = log[F(H j 9) /F(H/3)] 



(6) 



If an electron density of n c = 10 4 cm 3 is assumed, and T c = 
10 4 K, the Balmer-line intensity ratio, F(Ha)/F(H/3) = 2.85. 
Combining this factor with equations J5} & J6j gives an expres- 
sion for observed Ha extinction in terms of the observed Ha flux 
(mWm -2 ) and the (catalogue) radio flux (mJy): 



C a = log[2.855' y (6cm)/3.10 x 10 12 F(Ha) 



(7) 



In order to compare observed versus catalogue extinction val- 
ues, the latter can be derived from just the catalogue Ha and H/3 
flux values. We assume that the F(Ha)/_F(H/3) ratio in the cata- 
logue is accurate, even if the absolute fluxes are in doubt. We start 
with the full form of equation|6| 



Cp = log[F(H/9) /F(H/3)] = A iS 6iE B -v/2.5 



(8) 



where A\ is a wavelength dependent extinction coefficient. 
It follows that 



C a = log[F(Ha) /T(Ha)] = A 6561 E B - V /2.5 



(9) 



Combining equations|S|&|9|with the Balmer-line intensity ra- 
tio of 2.85 gives an expression for calculating Ha extinction values 
in terms of catalogue Ha and H/3 flux values: 



C aopt = Blog[T(Ha) cat /2.85F(H/3) cat ], 



(10) 



where B = ^4656i/(^4486i — ^656i) = 2.28 (using data from SP 
table V-l). 

An alternate method of calculating the wavelength 
depen dent constant B in equa tion 1101 is provided by 
ICardelli. Clayton & Mathisl j 19891 hereafter CCM) in 
terms of what they consider a more fundamental ex- 
tinction law A(\)/A(V), where A(X) is the abso- 
lute extinction at A and A(V) is the visual extinc- 
tion. Based on empirical data, CCM provide a mean 
Ry [= A(V)/E(B - V)] extinction law with wavelength 
dependent coefficients (x = A~ x /im _1 ) that takes the form 



(A(A)/A(V)) = a(x) + b(x)/Rv, 



(11) 



where Ry is the ratio of total to selective extinction. Dividing the 
expression for B by A(V) and substituting the expression a(x) + 
b(x)/Rv gives 



B 



{a a + b a /R\ 



(a/3 + bp/Ry) — (a a + b a /Rv) ' 



(12) 
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Figure 7. (a) Compai'ison of observed Ha extinction values, C a = 
log [2.85 S„ (6cm)/3.10 X 10 12 F(Ha)] and catalogued Ha extinc- 
tion values, C Qopt = B (Ry) log [F(Ha) cat /2.85 -F(H/3) cat ] where 
-B(3.1) = 2.36, for Bulge subset B objects. Diamonds denote objects 
observed on the first night, triangles the second, (b) Comparison of Ry 
(equation ll4l and observed Ha extinction values. 004.0— 11.1 (M3-29) and 
355.1-06.9 (M3-21) are excluded. 



where a a — a(l/AH a ), a/3 = (i(1/Ah/3), etc. Using Ry = 3.1 
and appropriate values provided by CCM (eqs. 3a & 3b) for this 
polynomial in x gives B = 2.36, which is slightly higher than SP's 
value of 2.28. However, B is now dependent on the value of Ry so 
equation fTol takes the Ry dependent form: 



C" Qopt = B(i*v)log[F(Ha) cat /2.85F(H/?)„t] 



(13) 



This enables Ry to be calculated in terms of the observed extinc- 
tion C a (equation[7) and catalogued Ha and H/3 flux ratios: 



Ry 



C a (b a -b ) + Db a 



C a (a,{3 - a a ) - Da a ' 
where D = log[F(Ha) ca t/2.85F(H/3) cat ]. 



(14) 



5.2 Results 

Our observed (equation and catalogue (equation 1 101 extinction 
values are listed in columns nine and ten of Table|5| At high bright- 
ness temperatures, the optical depth of PNe reduces observed ra- 
dio emission, so the observed extinction values are likely to be un- 
reliable, as are those calculated from radio fluxes below 10 mjy 
(indicated in column 13). Therefore Fig. Q(a) only compares ob- 
served and catalogue extinction values for a subset of Bulge ob- 
jects (hereafter subset B) with 100 mJy > S v > 10 mJy and 
T b < 10 3 K, where T b = 70.6255 GHz/6>i inc 2 (mJy/"). It can 
be seen that the values determined from our Ha flux and cat- 
alogued radio flux tend to be lower than those calculated with 
Ry =3.1 and catalogue Ha and H/3 fluxes. Subs et B conforms to 
the Bulge criteria of lBensbv & Lundst rom 1 2001 ) with the excep- 
tion of 004.0-11.1 (M3-29) where b = -11.1° and 350.9+04.4 
(H2-1) where Ssghz = 61 mJy. 

Calculated values for Ry (equation 1 141 are listed in column 
11 of Table Fig. 0(b) compares Ry with observed Ha extinc- 
tion values for subset B. The value of Ry is seen to be highly 
sensitive to differences between C a and C Q!opt with values > 10 
for objects 004.0-11.1 (M3-29) and 355.1-06.9 (M3-21) not in- 
cluded in the figure for the sake of clarity (regardless of method, 
calculated extinction values for these two objects are very low with 
AF(Ha) = 0.1, suggesting a possible error in the catalogue ra- 
dio values). Apart from the value for 009.4-09.8 (M3-32) and the 
above two objects, there appears to be a correlation between in- 
creasing values of Ry and observed Ha extinction. 

Because of the sensitivity of Ry to differences between ob- 
served and catalogue extinction values, and the uncertainty in these 
values, a reasonable error ARy has been calculated with C a ± 0.1 
using 



ARy = Ry(C a + 0.1) - Ry(C a - 0.1). 



(15) 



Using ARy as the error in Ry, the mean value (Ry) for subset B 
is found to be 2.0 (indicated on Fig.|7J. Values for ARy are listed 
in column 12 of Table [5] For subset B a comparison of Ry with 
catalogue 6 cm radio values shows no evidence of a correlation 
between low values of Ry and low radio fluxes. In fact, radio flux 
values would have to increase by a factor ~ 3 in order for (Ry) ~ 
3.1. 

Fig. [8] compares Ry with Galactic distribution for subset B 
(again excluding 004.0— 1 1.1 and 355.1—06.9). For the line of sight 
towards 351.1+04.8, the three objects close together on the sky 
(H2-1, Ml-19 and M2-5) have, within their uncertainties, (Ry) = 
1.2. It is worth noting that the two Sagittarius dwarf galaxy objects 
004.8-22.7 (He2-436) and 006.8-19.8 (Wray 16-423) have Ry = 
3.1 within their uncertainties. 

Objects 352.6+00.1 (Hl-12) and 352.8-00.2 (Hl-13) have 
not been included in subset B because of th eir very high r adio flux 
and being in the Galactic plane, leas well & Havneslfl987l) identify 
these positions with Hll regions with Ss ghz of 2.1 and 2.8 Jy re- 
spectively, giving higher values of Ry (2.4 and 3.5 respectively). 

Values for i?v C at and Ai?v C at were also calculated using 
catalogue Ha flux values in place of our observed values. The 
mean catalogue value (i?vcat) for subset B was found to be 2.2 
and for the line of sight towards 351.1+04.8 it was 1.4. A com- 
parison of observed and catalogue Ry values shows good agree- 
ment for subset B, apart from values for 009.4-09.8 (M3-32), 
352.1+05.1 (M2-8) and 003.6+03.1 (M2-14). This suggests that 
in addition to our newly observed Ha flux values, existing cat- 
alogue flux values also predict that interstellar extinction toward 
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Figure 8. Comparison of Ry with Galactic distribution for Bulge subset B 
objects. Relative Ry is indicated by dot size with Ry = 2.0 scale indicated 
at top right. 004.0-11.1 (M3-29) and 355.1-06.9 (M3-21) are excluded. 

the Bulge is lower than that corresponding to the standard extinc- 
tion curve Ry = 3.1. A comparison of observed and catalogue 
Su/Ha extinction values also shows good agreement for subset B, 
apart from values for 352.1+05.1 (M2-8) and 003.6+03.1 (M2- 
14), which have differences between observed and catalogued Ha 
flux values, AF(Ha) of 0.4 and 0.8 respectively. Excluding the 
above three objects from subset B has no effect of the calculated 
values of {Ry) and {Ry C at). There is no evidence for systematic 
differences due to time of observation or position in the sky for 
these three objects. 



6 DISCUSSION: EXTINCTION 

6.1 Evidence for anomalous extinction 

In an optical and UV analysis of Bulge PNe, 
IWalton. Barlow & Clesd 1 19931) found that Ry = 2.29 ± 0.50, 
concluding tha t the reddening law is steeper towards the Bulge. 
IStasinska etafl Il992t hereafter STAS) and iTvlenda et alJ <1992j 
hereafter TASK) discuss methods for determining the interstellar 
extinction of Galactic PNe and compare the two most frequently 
used; the Balmer decrement Ha/H/3 ratio method and the radio/H/3 
flux ratio method. Like all methods, the Balmer decrement method 
depends on the adopted extinction curve (or mean interstellar 
extinction law). In addition, it requires good quality optical 
data acquired by filter photometry, or spectrophotometry that 
encompasses the full extent of the nebula. 

In principle, the radio/H/3 flux ratio method is directly capa- 
ble of determining the actual value of H/3 extinction, Cp. How- 
ever, optical thickness of the nebula can lead to an underestima- 
tion of C/3. In addition, the predicted H/3 flux depends on adopted 
electron temperature and He ionization. Blending of Hell Pick- 
ering 8 emission with H/3 flux also adds a (very) small error. 
Interferometric radio data is preferred over single dish measure- 
ments beca use of the prob l em of confusion with other sources. 
IZiilstra. Pottasch & Bignelll Jl989l) have noted that their VLA 
fluxes are systematically lower than the results from single dish 



observations, due to possible cutting off of emission from fainter, 
extended regions of PNe. However, as noted previously, for sub- 
set B these 6 cm radio fluxes would have to be underestimated 
by a factor of ~ 3 for (Rv) ~ 3.1. Allowing that radio fluxes 
below 10 mjy are often underestimated for many faint PN (e.g. 
IPottasch & Ziilstrd Il994l) . strong radio sources still have C op t 
higher than C ra d. 

TASK compared the extinction values C ra d and C op t for 
Bulge and Disk PNe (their figs. 6 to 11), and showed that, for all 
PNe, Crad becomes systematically smaller than C op t for increasing 
extinction. They proposed that one of the possible reasons (as men- 
tioned by Cahn, Kaler & Stanahellini 1992) may lie in the value of 
Ry adopted in the extinction curve. They estimated that Ry — 2.7 
would bring C op t to a rough agreement with C ra d. 

KooDen & Vergelv ( 1998) analysed the optical extinction of 
27 1 Bulge PNe, and found (see their figs. 1 and 2) a mean extinction 
around C op t = 2, and a systematic decrease of the extinction (from 
3.5 to 0.5) with increasing Galactic latitude, as well as a genuine 
scatter about the relation. Both are accounted for by a model of 
small dust clouds randomly distributed in an exponential disc. They 
found an extinction Ay — 1.4 in the plane, and = 27 to the centre, 
in agreement with far IR-studies (COBE results). 

The standard extinction law for the diffuse interstellar medium 
is based on observations of O and B stars within a few kpc, but these 
stars must suffer sufficient extinction to allow measurement of their 
reddening. This selection effect favours denser regions of the in- 
terstellar medium in the Solar neighbourhood. For PNe the lines 
of sight are longer and could cross a more dilute medium while 
still giving measurable extinction. A lower value of Ry (< 2.5) 
could be due to a different size distribution of dust particles (i.e. 
smaller) in these more dilute regions. For the line of sight to- 
ward HP 210121 obscur ed by the high-latitude cloud DBB 80, 
lLarson et alJ ll996l l20od> find Ry = 2.1 ± 0.2, which they at- 
tribute to an abundance of small dust grains. 

More recently Udalski 1200 3l) presents results of analysis of 
photom etric data of Bulge fields from the OGLE II microlensing 
survey lUdalski etal1E 002). He shows Ry in general to be much 
smaller than that corresponding to the standard extinction curve 
Ry = 3.1, and that it varies considerably along different lines of 
sight. IPopowskl fcOOOl) has also suggested that interstellar extinc- 
tion toward the Bulge might be anomalous. Using CCM's extinc- 
tion law model, Udalski calculates values for Ry ranging from 1.75 
to 2.50 in four small regions of the Bulge. Assumed mean bright- 
ness, metallicity distribution, and large differences in the mean dis- 
tance of red clump giants are all dismissed as unlikely reasons for 
these highly variable values of Ry. 

Our results confirm these previous findings, and we find evi- 
dence for even lower vales of Ry. We believe our flux and extinc- 
tion determinations to be accurate because: our observed Ha fluxes 
are based on filter photometry; catalogue Ha and H/3 flux ratios are 
normally based on single spectra and are unlikely to suffer from 
systematic errors (even if the absolute fluxes are in doubt); and al- 
though catalogue radio fluxes are less certain, we have used a more 
reliable subset of Bulge objects (S„ > 10 mjy and T b < 10 3 K) for 
calculating values for Ry and ARy. However, radio fluxes should 
be remeasured to obtain the most accurate value for Ry. Our mean 
value of (Ry) = 2.0 provides further evidence that most lines 
of sight towards the Bulge cross a less dense medium containing 
smaller dust grains. It may be the case that the standard extinction 
law for the diffuse interstellar medium only applies to the local re- 
gion within a few kpc of our Sun. 
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6.2 Extinction and the warm ionized medium 

The PNe in our Bulge subset B sample have extinction up to Ay w 
5. The usual conversion between extinction and hydrogen column 
density, e.g., 

N(H) ^Ay2x 10 22 cm" 2 (16) 

gives column densities up to 10 23 cm -2 . For a line of sight of 8 
kpc, the average density becomes tih ^ 4cm -3 . 

Components of the interstellar medium are: Molecular clouds; 
Diffuse clouds; Warm neutral medium (WNM); Warm ionized 
medium (WIM); and Hot ionized medium (HIM). The warm neu- 
tral medium (WNM) has a typical density of 0.1-lOcm -3 , and 
scale height around 220 pc. The warm ionized medium (WIM) has 
a typical density of 0.3-10 cm" 3 in the Galactic plane, with a scale 
height of 1 kpc iBoul ares & Coxlll990l) . The HIM has density be- 
low 0.01 cm -3 and a scale height of 3 kpc. The WNM, WIM and 
HIM account for most of the volume in the Galactic disc. The 
WNM is absent near the Sun (within lOOpc in most directions): 
the local volume is filled with a HIM with embedded cold clouds. 

Our extinction values of A\r ^ 5 mag suggest that the lines 
of sight do not intercept dense molecular clouds. The density in the 
HIM is too low to give appreciable extinction even out to the Bulge. 
The extinction to the Bulge PNe is likely to arise from a mixture of 
WIM and WNM. The line of sight at higher latitudes leaves the 
WNM within 2-3 kpc (no Bulge PNe are known within 2 degrees 
of the plane): the remainder will be mainly within the WIM. 

The relation between Ry and C a depicted in Fig. 0(b) sug- 
gests that the lowest R values correspond to the lower density ISM 
component. A large fraction of the lines of sight fall in WIM- 
dominated regions above/below the disc. The WNM component is 
more likely to cause the variable amount of additional extinction. 
The fact that three PNe located within a degree of each other show 
very similar - and very low - Ry shows that the low-density struc- 
tures are similar for line of sights within approximately 100 pc: this 
provides support for the identification of this component with the 
WIM. If correct, the distribution of Ry in Fig. 0(b) suggests that 
the WIM has Ry » 1.5, and the WNM Ry « 2.5. The latter value 
is close to the lowest local values found by CCM. 

The anomalous extinction towards the Bulge therefore suggest 
that the WIM shows smaller dust grains than found in the WNM or 
molecular clouds. Possible explanations include dust-grain evolu- 
tion, and the effects of supernovae. The ISM cycles between high- 
density and low-density phases. If the dust is slowly destroyed dur- 
ing the low-density phases, the size of the dust grains may measure 
the time since the last high density phase. The large scale height of 
the WIM suggests its contents may cycle much slower. 

Supernova-driven shocks can also destroy dust grains. Super- 
novae are more dominant in the inner galaxy <Heileslll987l) : The 
low-density medium in the inner Galaxy could therefore show sys- 
tematically smaller dust grains. This would limit the applicability 
of the standard Galactic extinction curve to quiescent regions of the 
Galaxy, with steeper curves found in energetic environments. 

6.3 Grain sizes 

Although different from the locally derived values of Ry, it is pos- 
sible to show that our low values derived above can be reproduced 
using dust models. 

The Gr eenberg unified dust model <Greenburg&Tlll996t 
iLi & Greenberg! \l 9971) is capable of generating a variety of val- 
ues of Ry, which encompass the observational value derived in 



the present work, by varying the turnover radius, a key parameter 
of the mod el. 

lLi&Greenberdi!99l) use a three component model to com- 
pute the entire interstellar extinction law: large composite grains, 
small carbonaceous grains and a population of PAH molecules. 
In the present work, we consider only the large composite grains, 
which consist of a silicate core and a refractory organic mantle. Ad- 
dition of the smaller components would tend to increase ultraviolet 
extinction relative to that in the visible, and thus drive the model 
results in the present work to smaller values of Ry. The size dis- 
tribut ion function used for the composite grains is iLi & Greenberg! 
fl99l) 

n(a) = n exp { ~ 5 ( ° fl " C ) } ( 17 ) 

where n(a) is the number density of grains of total radius a. 
These composite grains have a silicate core of radius a c , and a 
turnover radius, a m . The exponent, q, ensures that the distribu- 
tion of grains falls rapidly for sizes la rger than q m . In the p resent 
work we use values of q = 2, as in lLi & Greenberg! j 1997l> . and 
set a c = 0.042 /im which was the radius of infinite cylinders in 
Li & Greenberg, but spheres in the present work. The value of a m 
was varied to generate a tabl e of theoreti cal values of Ry . The core 
material used by Greenburg & Li 1 1996) was based on values for 
an amorphous olivine l Dorschne r et alll995t) . In the present work 
we replace t he olivine by am orphous forsterite, defined by optical 
constants (in Scott & Dulev 1996), on the grounds that the olivine 
contains iron, whilst interstellar silicate is now believ ed to be domi- 
nated by magnesium-rich minerals (see, for example, Kim ura et alJ 
2003). Opti cal constants for the organic refractory mantle were 
taken from lGree nburg T3C996). 

Fig. |5| (left panel) shows the optical scattering and absorp- 
tion coefficients, and the extinction law, for the adopted 'standard' 
model, with one particular value of the turnover radius, a m — 
0.08 fim. In Fig. [5] (right panel), we plot the value of _Rv as a 
function of the turnover radius, a m . The theoretical values shown 
encompass all values of Ry from the smallest-but-one of the ob- 
servational figures obtained in the present work, through the rest 
of the present samples and the more traditional values of ~ 3.1 
expected for dust in diffuse clouds, to values in excess of 5, ex- 
pected for dense, cold, clouds. We note that the key parameter is 
the amount of grain material found in large grains (with radii larger 
than ~ 0.1 fjm). 

We also note that, for any value of a m , a smaller value of Ry 
may be obtained by adding the populations of small carbonaceous 
grains and PAHs, whilst a larger value is computable by overlaying 
the composite grains with layers of ices, as is likely in cold clou ds. 
For a comparison of ice-coated and bare grains, see lGravH200ll) . 



7 CONCLUSIONS 
7.1 Angular diameters 

We have calculated angular diameters for 69 PNe by deconvolv- 
ing the point spread function of field stars from Gaussian FWHM 
profiles. A correction factor, 7, was applied to determine the true 
diameter. 7 was approximated by a simple analytic function and 
diameters were calculated using two models: a constant emissiv- 
ity shell 0.8, #sh e il, ar, d photoionization line emission, #ii ne . We 
present the mean value derived from the two models, # mea n, as 
our best estimate. For some objects the low intensity contour plots 
revealed an elliptical structure that was not always apparent from 
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Figure 9. Left: A plot, between 0.3 and 1 .0 microns, of the absorption, scattering and extinction cross-sections for a Greenberg model with a turnover radius of 
0.08 microns, q = 2, and core radius equal to 0.042 microns. The cross-sections are measured relative to an area of 2.147 X 10 -16 m~ 2 ; Right: Calculated 
values for Rv for a Greenberg dust model, as function of the turnover radius for the size distribution. 



the FWHM measurements. Axis averaged mean angular diameters 
of fully resolved objects ranged from 2.8 ± 0.1 to 12.7 ± 1.2 a rc 
seconds. 18 Ha images were restored with the Richards orllLucvl al- 
gorithm and restored diameters tended to be around one arc second 
smaller below 4 arc seconds. 



7.2 Fluxes 

The total flux detected through the filters was calculated for 70 PNe. 
Doppler shifted Ha, [O III] (and [N II]) wavelengths observed at the 
telescope were computed from radial velocities, and filter transmis- 
sion coefficients applied to arrive at actual flux densities. For Ha 
and [O III] there is little evidence of any systematic differences be- 
tween observed and catalogued flux values. The accuracy of filter 
response data was considered the most significant source of error 
and estimated at between ±5 to ±10 per cent. 

7.3 Extinction 

Observed Ha extinction in the direction of the Bulge was deter- 
mined using the radio to Hf3 flux ratio, the expected to observed 
H/3 flux ratio and the Balmer-line intensity ratio, giving an ex- 
pression in terms of observed Ha flux and catalogue radio flux. 
Catalogue Ha extinction values were also derived from catalogued 
Ha and H/3 flux values in terms of the fundamental extinction law 
A(X) /A( V) by means of an Rv dependent function. We compared 
observed and catalogue extinction values for a subset of Bulge ob- 
jects and found that values determined from observed Ha flux and 
catalogued radio flux tended to be lower than those calculated with 
Rv — 3.1 and catalogue Ha and H/3 fluxes. A method for deter- 
mining i?v was derived in terms of observed extinction and cata- 
logued H/3 and Ha flux values. We found Rv to be highly sensitive 
to differences between observed and catalogue extinction values. 
There also appears to be a correlation between increasing values of 
both Rv and observed Ha extinction, which is in agreement with 
the findings of STAS discussed above. Estimating a reasonable er- 
ror in Ry, we find observed {Rv} = 2.0, which fits within the 



range calculated by Udalski and agrees with his conclusion that to- 
ward the Bulge interstellar extinction is steeper than Rv = 3.1. 
Using catalogue Ha and S v flux values we find good agreement, 
with (Rvcat) = 2.2. For one small region in our sample we find 
(Rv) as low as 1.2, which suggests that Rv does indeed vary con- 
siderably along different lines of sight. We have shown that our 
values of Rv can be reproduced using dust models with a turnover 
radius of 0.08 microns. A lower value of Rv may also be obtained 
by adding the populations of small carbonaceous grains and PAHs. 
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